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METHODS FOR GEOMECHANICAL FRACTURE MODELING 
Background of the Invention 

[0001] The present invention relates generally to methods for designing and optimizing 
the number, placement, and size of fractures in a subterranean formation and more particularly to 
methods that account for stress interference from other fractures when designing and optimizing 
the number, placement, and size of fractures in the subterranean formation. 

[0002] One method typically used to increase the effective drainage area of well bores 
penetrating geologic formations is fracture stimulation. Fracture stimulation comprises the 
intentional fracturing of the subterranean formation by pumping a fracturing fluid into a well bore 
and against a selected surface of a subterranean formation intersected by the well bore. The 
fracturing fluid is pumped at a pressure sufficient that the earthen material in the subterranean 
formation breaks or separates to initiate a fracture in the formation. 

[0003] Fracture stimulation can be used in both vertical and horizontal wells. Fracturing 
horizontal wells may be undertaken in several situations, including situations where the formation 
has: 

1 . restricted vertical flow caused by low vertical permeability or the presence of 
shale streaks; 

2. low productivity due to low formation permeability; 

3. natural fractures in a direction different from that of induced fractures, thus 
induced fractures have a high chance of intercepting the natural fractures; or 

4. low stress contrast between the pay zone and the surrounding layers. 

In the fourth case, a large fracturing treatment of a vertical well would not be an acceptable 
option since the fracture would grow in height as well as length. Drilling a horizontal well and 
creating either several transverse or longitudinal fractures may allow rapid depletion of the 
reservoir through one or more fractures. 

[0004] Shown in Figure 1 is an example of a well bore, represented generally by the 
numeral 100, comprising a generally vertical portion 102 and two laterals 104 and 106. The 
generally vertical portion 102 is drilled in a generally vertical direction, and the laterals 104 and 
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106 are disposed at angles 108 and 110, respectively to the vertical portion 102. The well bore 
100 is referred to as a horizontal well because it has one or more laterals (in the case of well 100, 
laterals 104 and 106). Typically, only the laterals 104 and 108 are open for production in a 
horizontal welL If the well 100 only had a generally vertical portion 102, it would be referred to 
as a vertical well. Typically, all production in a vertical well comes from the generally vertical 
portion 102. 

[0005] Shown in Figure 2 is a perspective view of the well bore 100 comprising lateral 
104. The lateral 104 comprises three fractures 202, 204 and 206. Depending on the orientation of 
the lateral 204 to the direction of minimal stress, the fractures 202, 204 and 206 may be transverse 
or axial fractures. If the lateral 104 is drilled in direction of minimal stress, then the fractures 202, 
204 and 206 are orientated perpendicular to the direction of minimal stress and are referred to as 
transverse fractures. If the lateral 104 is drilled perpendicular to the direction of minimal stress, 
then the fractures 202, 204 and 206 are orientated parallel to the direction of minimal stress and 
are referred to as axial fractures. 

[0006] Each of the fractures 202, 204 and 206 typically has a narrow opening that 
extends laterally from the well bore. To prevent such opening from closing completely when the 
fracturing pressure is relieved, the fracturing fluid typically carries a granular or particulate 
material, referred to as "proppant," into the opening of the fracture and deep into the fracture. 
This material remains in each of the fractures 202, 204 and 206 after the fracturing process is 
finished. Ideally, the proppant in each of the fractures 202, 204 and 206 holds apart the separated 
earthen walls of the formation to keep the fracture open and to provide flow paths through which 
hydrocarbons from the formation can flow into the well bore at increased rates relative to the flow 
rates through the xmfractured formation. Fracturing processes are intended to enhance 
hydrocarbon production from the fractured formation. In some circumstances, however, the 
fracturing process may terminate prematurely, for a variety of reasons. For example, the "pad" 
portion of the fracturing fluid, which is intended to advance ahead of the proppant as the fracture 
progresses, may undesirably completely "leak ofF' into the formation, which may cause the 
proppant to reach the fracture tip and create an undesirable "screenout" condition. Thus, properly 
predicting fracture behavior is a very important aspect of the fracturing process. 
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[0007] In the past, fracturing typically took place in well bores that were cased and 
perforated. The total number of fractures was a limited number per lateral in the case of 
fracturing horizontal wells and the fractures had sufficient space between each other such that 
stress interference between the fractures was minimal With the advent of new fracturing 
technologies such as SURGIFRAC provided by Halliburton Energy Services, fractures may be 
placed in open hole well bores. Furthermore, it is now feasible and cost-effective to place many 
more fractures in a well bore. When many fractures are induced in a well bore, the 
geomechanical stress caused by fractures on each other can no longer be ignored. Current 
fracturing modeling methods, however, do not accoxmt for geomechanical stresses caused by one 
fracture on another. 
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Summary of the Invention 

[0008] The present invention relates generally to methods for designing and optimizing 
the number, placement, and size of fractures in a subterranean formation and more particularly to 
methods that account for stress interference from other fractures when designing and optimizing 
the number, placement, and size of fractures in the subterranean formation. 

[0009] One embodiment of the present invention includes a method of optimizing a 
number, placement and size of fractures in a subterranean formation, including the steps of: 
determining one or more geomechanical stresses induced by each fracture based on the 
dimensions and location of each fracture; determining a geomechanical maximum number of 
fractures based on the geomechanical stresses induced by each of the fractures; and determining a 
predicted stress field based on the geomechanical stresses induced by each fracture. 

[0010] Another embodiment of the present invention includes a computer implemented 
method for optimizing a number, placement and size of fractures in a subterranean formation, 
including the steps of: determining one or more geomechanical stresses induced by each fracture 
based on the dimensions and location of each fracture; determining a geomechanical maximum 
number of fractures based on the geomechanical stresses induced by each of the fractures; and 
determining a predicted stress field based on the geomechanical stresses induced by each fracture. 

[0011] Another embodiment of the present invention includes a method of fracturing a 
subterranean formation, including the step of optimizing a number, placement and size of 
fractures in the subterranean formation, the step of optimizing including the steps of: determining 
one or more geomechanical stresses induced by each fracture based on the dimensions and 
location of each fracture; determining a geomechanical maximum number of fractures based on 
the geomechanical stresses induced by each of the fractures; determining a predicted stress field 
based on the geomechanical stresses induced by each fracture. 

[0012] The features and advantage of the present invention will be readily apparent to 
those skilled in the art upon a reading of the description of the preferred embodiments which 
follows. 



4 



Brief Description of the Drawings 

[0013] The present invention is better understood by reading the following description of 
non-limitative embodiments with reference to the attached drawings wherein like parts of each of 
the several figures are identified by the same referenced characters, and which are briefly 
described as follows: 

[0014] Figure 1 depicts a side cross-sectional view of a subterranean well bore having 

laterals. 

[0015] Figure 2 is a perspective view of a well bore with a lateral having fi-actures. 

[0016] Figure 3 illustrates a process flow diagram fi-om an exemplary method of the 
present invention for creating a fi'acture layout. 

[0017] Figure 4 illustrates a process flow diagram fi-om an exemplary method of the 
present invention for estimating a cost-effective number of fi'actures. 

[0018] Figure 5 illustrates a process flow diagram from an exemplary method of the 
present invention for estimating a geomechanical maximum number of fractures. 

[0019] Figure 6 illustrates a process flow diagram from an exemplary method of the 
present invention for modeling a fracture. 

[0020] Figure 7 is a graphical representation of the principal components of stress 
induced by a semi-infinite fracture versus dimensionless distance. 

[0021] Figure 8 is a graphical representation of the principal components of stress 
induced by a permy-shaped fracture versus dimensionless distance. 

[0022] Figure 9 is a graphical representation of the principal components of stress 
induced by a semi-infinite fracture and a penny-shaped fracture versus dimensionless distance. 

[0023] Figure 10 illustrates a process flow diagram from an exemplary method of the 
present invention for creating a fracture layout. 

[0024] Figure 1 1 illustrates a process flow diagram from an exemplary method of the 
present invention for determining whether modeled fractxu-es fail. 

[0025] Figure 12 depicts a side cross-sectional view of a subterranean well bore wherein 
fluid may be injected, and the results of such injection monitored, according to an exemplary 
embodiment of the present invention. 
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[0026] Figure 13 illustrates a process flow diagram from an exemplary method of the 
present invention for fracturing based on a fracture layout. 

[0027] Figure 14 is a graphical representation of instantaneous shut-in pressure versus 
time for the creation of six fractures. 

[0028] Figure 15 is a graphical representation of instantaneous shut-in pressure and true 
vertical depth versus measured depth for the creation of twenty fractures. 

[0029] Figure 16 is a graphical representation of instantaneous shut-in pressure and true 
vertical depth versus measured depth for the creation of twelve fractures. 

[0030] It is to be noted, however, that the appended drawings illustrate only typical 
embodiments of this invention and are therefore not to be considered Umiting of its scope, as the 
invention may admit to other equally effective embodiments. 
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Detailed Description Of The Invention 

[0031] The present invention relates generally to methods for designing and optimizing 
the number, placement, and size of fractures in a subterranean formation and more particularly to 
methods that account for stress interference from other fractures when designing and optimizing 
the number, placement, and size of fractures in the subterranean formation. The present invention 
may be applied to vertical or horizontal wells. Furthermore, the present invention may be used on 
cased well bores or open holes. 

[0032] Figure 3 depicts a flow chart of an exemplary embodiment of the methods 
according to the present invention. In step 302, the method determines a cost-effective number of 
fractures. In step 304, the method determines a geomechanical maximum number of fractxxres. In 
step 306, the method determines whether the cost-effective number of fractures or the 
geomechanical number of fractures is limiting. If the cost-effective number of fractures is 
limiting (e.g,, if the method determines that geomechanically the formation can sustain more 
fractures than are cost-effective) then the method proceeds to step 308 where it creates a fracture 
layout based on the cost-effective number of fractures. If the geomechanical maximum number of 
fractures is limiting (e.g., if the method determines that geomechanically the formation can 
sustain less fractures than are cost-effective) then the method proceeds to step 310 where it 
creates a fracture layout based on the geomechanical maximum munber of fractures. 

[0033] Referring now to Figure 4, step 302 of Figure 3, in which the method according to 
the present invention determines the cost-effective number of fractures, is shown in greater detail 
In step 402, the method sets the cost-effective number of fractures to zero. In step 404, the 
method estimates an increase in production for a next modeled fracture. Referring to Figure 2, 
each of the fractures 202, 204, and 206 has an associated increase in production. Typically, the 
associated increase in production of a next modeled fracture is smaller than the increase in 
production associated with a previously modeled fracture. The increase in production of each 
additional fracture may be calculated based on any conventional method. In an exemplary 
embodiment of the present invention the method may consider some or all of the following 
criteria to determine the increase in production for the next fracture: physical properties of the 
formation (e.g., horizontal and vertical permeability, whether anisotropy is present, whether the 
formation if homogenous or heterogeneous, vertical lithological definitions including layers and 
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shale streaks, and a leak of coefficient), physical properties of the reservoir {e,g., pressure, 
porosity, height, temperature, formation compressibility, fluid saturation, a type of fluid in the 
reservoir, and properties of the fluid in the reservoir), a definition of the stress field {e.g., a 
minimum horizontal stress in a pay zone and surrounding zones and a stress orientation of the 
formation), and mechanical properties of the rock in the formation (e,g,, a Young's modulus due 
to the rock and a Poisson's ratio due to the rock). 

[0034] Returning to Figure 4, in step 406, the method estimates the cost of the next 
modeled fi-acture. The cost of each additional fi-acture is determined by adding all costs 
associated with the next modeled fi-acture. hi step 408, the method calculates the cost-benefit 
ratio by dividing the estimated cost associated with the next modeled firacture by the estimated 
increase in production associated with the next modeled fracture. Li block 410, the method 
determines if the cost-benefit ratio for the next modeled firacture is greater than a maximmn cost- 
benefit ratio. The maximum cost-benefit ratio may be set by the user on a case-by-case basis or 
may be a default value. If the cost-benefit ratio for the next modeled fracture is greater than the 
maximum cost-benefit ratio the method proceeds to step 412, where the cost-effective number of 
fractures is retumed. If the cost-benefit ratio for the next modeled fracture is not greater than the 
maximum cost-benefit ratio then the method proceeds to step 412, where the cost-effective 
number of fractures is increment by one and the routine is repeated by returning to step 404. 

[0035] The methods of the present invention may use metrics other than cost-benefit ratio 
for optimizing the number of fractures. For example the method of the present invention may use 
other financial parameters including a net present value (NPV) of each fracture, a pay-out time of 
each of the fractures, or other financial parameters of creating each of the fractures. 

[0036] Referring now to Figure 5, step 302 of Figure 3, in which the method according to 
the present invention determines the geomechanical maximum number of fractures, is shown in 
greater detail. In step 502, the method sets the geological maximum number of fractures to zero. 
In step 504, the method determines an initial stress field of the well bore in the geological 
formation. In step 506, the method determines if a next modeled fracture will fail. If the next 
modeled fracture fails the method proceeds to step 508, where it retums the geological maximum 
number of fractures. If the next modeled fracture does not fail the method proceeds to step 510, 
where it models the next modeled fracture. In step 512, the method determines the new stress 
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field due to the placement of the next modeled fi-acture. In step 514, the method increments the 
geological maximum number of fi-actures by one and returns to step 506. 

[00371 In step 504, the method determines an initial stress field of the well bore in the 
geological formation. Referring to Figure 1, the initial stress field on well bore 100 may be input 
by the user or deteraiined by any conventional method including sampled data fi-om the formation 
including microfi-acturing test data, minifi-acturing test data, leak-off test (LOT) data, or logging 
data. In an exemplary embodiment of the present invention wavelet analysis is used to determine 
the stresses fi-om microfi-acturing or minifracturing test data. The method then determines the 
orientation of the vertical portion 102 to the initial stress field. The orientation of the vertical 
portion 102 may be input by the user or the method may determine the orientation of the vertical 
portion 102. In an exemplary embodiment of the present invention, the method determines the 
orientation of the well bore 102 by assuming that the well bore 102 will be placed parallel to the 
direction of maximum stress (overburden stress) in the initial stress field. If the method is 
determining the placement of fi-actures in a horizontal well, the method determines the orientation 
of one of the laterals 104 or 106 to the initial stress field. The orientation of the one of the laterals 
104 or 106 may be input by the user or may be determined by the method. In an exemplary 
embodiment of the present invention, the method determines the orientation of the one of the 
laterals 104 or 106 by assuming that the one of the laterals 104 or 106 will be orientated parallel 
to the direction of minimimi stress in the initial geological formation. 

[0038] Referring again to Figure 5, in step 506, the method of the present invention 
determines if the next modeled fracture will fail. The next modeled fracture will fail when it 
propagates in a tortuous patti, leading to higher fracture pressure and possibly to sand-out. For 
example, if a transverse fracture is placed in a lateral of a horizontal well bore, it will fail if it 
"turns" and begins to propagate in an axial direction. In another example, if an axial fracture is 
placed in a vertical well bore, it will fail if it "tums" and begins to propagate in a transverse 
direction. To predict if a fracture will fail, the method of the present invention calculates the 
geomechanical stresses at the point where the modeled fracture is initiated. To determine the 
point where the next modeled fracture will be initiated the method may receive input from the 
user or the method may determine the point where the next modeled fracture will be initiated 
automatically. In an exemplary embodiment of the present invention, the method assumes that 
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the modeled fractures are equidistant from each other. The method calculates the geomechanical 
stresses at the point where the next modeled fracture is initiated by summing the initial stress field 
and the stress fields caused by any previous modeled fractures. After this summation, the method 
determines which principal component of geomechanical stress is smallest at the point where the 
modeled fracture is initiated. Li the case of a transverse fracture in a lateral of a horizontal well 
bore, if the minimum stress is the vertical stress then the fracture is deemed to fail. In the case of 
an axial fracture in a vertical well bore, if the minimum stress is the horizontal stress the fracture 
is deemed to fail. 

[0039] Referring now to Figure 6, step 510, in which the exemplary method of the 
current invention models the next modeled fracture, is shown in greater detail. Li step 502, the 
method selects a model to use to model the fracture. The selection of one of the models may be 
accomplished with or without user intervention. In an exemplary embodiment of the present 
invention, the user manually selects a model to use for modeling the next modeled fracture and 
inputs the dimension of the fracture. In another embodiment of the present invention, there is a 
defauh fracture model used to model the next modeled fracture. In yet another embodiment of the 
present invention, the method will determine which model is most appropriate for modeUng the 
next modeled fracture based on the input characteristics of the next modeled fracture and 
previously modeled fractures (e.g., the distance between fractures, the size of the fracture, and the 
shape of the fracture). If the method chooses to model the next modeled fracture as a semi- 
infinite crack the method proceeds to step 604. If the method chooses to model the next modeled 
fracture as a penny-shaped fracture the method proceeds to step 606. If the method chooses to 
model the next modeled fracture using another geomechanical model the method proceeds to step 
608. 

[0040] Regardless of the method used to model the next modeled fracture, the method of 
the present invention may consider properties of the geological formation (e,g. type of material 
and presence of naturally occurring fractures) while modeling the next modeled fracture. In an 
exemplary embodiment of the present invention the method considers the presence of naturally 
occurring fractures in the geological formation. The presence of these fractures may reduce the 
stress induced by the previously modeled fractures on the next modeled fracture. 
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[0041] When modeling the next modeled fracture as a semi-infinite crack in step 604, the 
method of the present invention assumes that next modeled fracture is rectangular, with an infinite 
length, a finite height, and a width that is extremely small compared with the height and the 
length of the fracture. The height of the next modeled fracture may be input by the user or may be 
determined by the method. In an exemplary embodiment of the present invention, the method 
assumes that the modeled fractures have equal dimensions, and optimizes the size of the fractures 
to maximize the geological maximum number of fractures. Using these assumptions the method 
of the present invention calculates the stress field caused by the next modeled fracture using the 
following equations: 
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where: Ox, Oy, and Oz are the components of stress in the x, y, and z directions respectively; is 
the shearing stress; po is the internal pressure at the point where the fracture is initiated; H is the 

height of the fracture ;jL/ is Poisson's ratio; and where z = re'^ , z-—H = , z-\-—H = r^e'^' . 

The method also records a predicted fracturing pressure associated with the next modeled 
fracture. In an exemplary embodiment of the present invention, the predicted fracturing pressure 
is equal to the internal pressure. 

[0042] Referring now to Figure 7, depicted is a graphical representation of the change in 
the three components of the principal stresses (Ox, Cy, and oz) versus the ratio L/H where L is a 
distance from the fracture along a line of symmetry and H is the height of the fracture. The line of 
symmetry is used because it represents the horizontal direction in case of creation of multiple 
fractures from a horizontal well. With respect to the coordinates of the fimctions plotted in Figure 
8, the x-direction is the direction perpendicular to the created fracture, the y-direction is the 
horizontal direction parallel to the fracture, and the z-direction is the vertical direction. 
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[0043] Referring again to Figure 6, when modeling the next modeled fracture as a penny- 
shaped fracture in step 606, the method of the present invention assumes that the next modeled 
fracture is circular shaped and has finite dimensions. The height of the next modeled fracture may 
be input by the user or may be determined by the method. In an exemplary embodiment of the 
present invention, the method assumes that the modeled fractures have equal dimensions, and 
optimizes the size of the fractures to maximize the geological maximmn number of fractures. 
Using these assumptions the method of the present invention calculates the stress field caused by 
the next modeled fracture using the following equations: 
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where: Or, Gz, and oe are the polar components of stress; r^^ is the shearing stress; po is the 
internal at the point where the fracture is initiated; z = re'^ , z-c = r^e'^ , and z + c = r^e'^^ , 
where the fracture extends from z=c to z=-c; and where a two-dimensional projection of the 
fracture is defined by the fimction 77^ = -b^ , where the origin of the coordinates is the edge of 
the fracture, ^ is the axis along the fracture, r| is the axis perpendicular to the fracture, 

= 5 cos ^ , and r| = 5 sin 4^ . The equations are provided in this coordinate set for brevity. One 
of ordinary skill in the art with the benefit of this disclosure can convert the coordinates and 
solve for Ox, Cy, and Oz. The method also records a predicted fracturing pressure associated with 
the next modeled fracture. In an exemplary embodiment of the present invention, the predicted 
fracturing pressure is equal to the internal pressure. 

[0044] Referring now to Figure 8, depicted is a graphical representation of the change in 
the three principal stresses (Ox, Oy, and Oz) versus the dimensionless distance L/H where L is the 
distance from the fracture and H is the diameter of the fracture for the penny-shaped fracture. 
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With respect to the coordinates of the functions plotted in Figure 9, the x-direction is the direction 
perpendicular to the created fracture, the y-direction is the horizontal direction parallel to the 
fracture, and the z-direction is the vertical direction. 

[0045] Referring now to Figure 9, depicted is a graphical representation of the change in 
minimum horizontal stress (the stress component perpendicular to the fracture) due to the creation 
of a semi-infinite fracture versus dimensionless distance from the fracture and the change in 
minimum horizontal stress due to the creation of a penny-shaped fracture versus dimensionless 
distance from the fracture. The dimensionless distance from the fracture is the ratio of the 
distance from the fracture versus the height or diameter of the fracture. 

[0046] Referring again to Figure 6, in step 608, the method according to the present 
invention may use other geomechanical models to model the next modeled fracture. In one 
exemplary embodiment of the present invention, the method may model the fractures as both a 
penny-shaped fracture (as in step 604) and as a semi-infinite fracture (as in step 602) and 
interpolate between the modeled stress fields (penny-shaped and semi-infinite) based on one or 
more properties of the next modeled fracture (e.g. the length of the next modeled fracture or the 
shape of the next modeled fracture) to determine a stress field for the modeled fracture. In an 
exemplary embodiment of the present invention the dimensions of the next modeled fracture are 
input by the user. In another exemplary embodiment of the present invention, the method 
assumes that the modeled fractures have equal dimensions, and optimizes the size of the fractures 
to maximize the geological maximum number of fractures. The method may assign a weight to 
the length and diameter/height of the fracture. In that case, stress field induced by a longer 
fracture will more closely resemble the stress field induced by a semi-infinite fracture than a 
shorter fracture, assuming all other dimensions of the longer and shorter fractures are equivalent. 
The method also records a predicted fracturing pressure associated with the next modeled 
fracture. In an exemplary embodiment of the present invention, the predicted fracturing pressure 
is equal to the internal pressure. 

[0047] Referring again to Figure 5, the method determines the new stress field due to the 
next modeled fracture. The method sums the initial stress field, the stress fields caused by 
previously modeled fractures, and the stress field case by the next modeled fracture. In an 
exemplary embodiment of the present invention, it is assumed that the medium is linearly elastic 



13 



and that the governing model of the stress field (comprising the differential equations, boundary 
conditions, and initial conditions) is linear, the principle of superposition is applicable. Thus, the 
method of the present invention may calculate the new stress field by summing the stresses caused 
by each of the fi-actures on the specific point in the formation. 

[0048] In another exemplary embodiment of the present invention, the method may 
calculate the stress field by using superposition and by adding the initial stress field, the stress 
fields caused by each of previously modeled fi"actures, and the next modeled fi-acture, 
sequentially. This has the effect of predicting a greater change in the minimum stress because 
each modeled fi-acture will be created against a higher minimum stress (due to the presence of the 
previously modeled stress fields). Because the minimum stress will be higher for each subsequent 
fi^acture, the internal pressure at the point where the subsequent fracture is initiated will be higher. 
Consequently, a higher fracturing pressure will be required to create each subsequent fracture 
overcome the internal pressure of the formation. The increase in will, in turn, lead to a greater 
change in the minimum stress caused by the next modeled fracture. 

[0049] The method may also calculate the new stress field due to the creation of fractures 
in multiple laterals of a single well. For example, referring to Figure 1, fractures may be initiated 
in laterals 104 and 106. The method may calculate the new stress field for fractures initiated in 
lateral 106 including the stress field induced by fractures 202, 204, and 206 (shown in Figure 2) in 
lateral 104. The method may also calculate the stress field due to adjacent well bores or fractures 
in adjacent well bores around well bore 102. 

[0050] Referring again to Figure 3, in step 308, the method uses the fracture layout based 
on the cost-effective number of fractures. The method creates the fracture layout such that it has 
the cost-effective number of fractures. The method may use any conventional method to produce 
the fracture layout. The fracture layout may be generated on a computer and output to a display 
device or printer. The fracture layout may be controlled by the input of the user or the method 
may determine the fracture layout automatically. In an exemplary embodiment of the present 
invention, the method will create the fracture layout so that the fractures are spaced equally from 
each other. The size of the fractures may be input by the user or the method may determine the 
size of the fractures automatically. 
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[0051] In step 310, the method uses the fracture layout based on the geomechanical 
maximum number of fractures. The method creates the fracture layout such that it has the 
geomechanical maximum number of fractures. The method may use any conventional method to 
produce the fracture layout. The fracture layout may be generated on a computer and output to a 
display device or printer. The fracture layout may be controlled by the input of the user or the 
method may determine the fracture layout automatically. In an exemplary embodiment of the 
present invention, the method will create the fracture layout so that the fractures are spaced 
equally from each other. The size of the fractures may be input by the user or the method may 
determine the size of the fractures automatically. 

[0052] Figure 10 depicts a flow chart of another exemplary embodiment of the methods 
of present invention. In step 302 (discussed above with respect to Figure 3) the method 
determines a cost-effective number of fractures. In step 1004, the method determines if the cost- 
effective number of fractures fail. If the cost-effective number of fractures fails, the method 
proceeds to step 1 104, where the cost-effective number of fractures is decremented by one and the 
method proceeds to step 1002. If the cost-effective number of fractures do not fail the method 
proceeds to block 308 (described above with respect to Figure 3) where it creates a fractvu-e layout 
based on the cost-effective number of fractures. 

[0053] Referring now to Figure 11, step 1002 of Figure 10, in which the method 
according to the present invention determines whether the cost-effective number fractures will 
fail, is shown in greater detail. In step 1 102, the method sets a geologically modeled number of 
fractures to zero. In step 504, the method determines the initial stress field of the well bore in the 
geological formation, as described with respect to Figure 5. In step 506, the method determines if 
the next modeled fracture fails, as described with respect to Figure 5. If the next modeled fracture 
does fail, the method proceeds to step 1 104 where it returns "Yes." If the next modeled fracture 
does not fail, the method proceeds to step 1 106. In step 1 106, the method determines if the 
geologically modeled number of fractures is equal to the cost-effective number of fractures. If the 
geologically modeled number of fractures is equal to the cost-effective nimiber of fractures the 
method proceeds to step 1 108, where it retums "No." If the geologically modeled number of 
fractures is not equal to the cost-effective number of fractures, the method proceeds to step 510 
where it models the next modeled fracture, as described with respect to Figure 5. In step 512, the 
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method determines the new stress field due to the next modeled fi*actiire, as described with respect 
to Figure 5. In step 1110, the method increments the geologically modeled number of fractures by 
one and returns to step 506. 

[0054] Figure 12 depicts a schematic representation of a subterranean well bore 1212 
through which a fluid may be injected into a region of the subterranean formation sxirrounding 
well bore 1212 such that physical property data (e.g., pressure signals, temperature signals, and 
the like) are generated. The fluid may be of any composition suitable for the particular injection 
operation to be performed. For example, where the methods of the present invention are used in 
accordance with a fracture stimulation treatment, a fracturing fluid may be injected into a 
subterranean formation such that a fracture is created or extended in a region of the formation 
surrounding well bore 1212 and generates pressure signals. The fluid may be injected by injection 
device 1201 (e.g., a pump). Physical property data such as pressure signals may be generated 
during subterranean injection processes, for reasons including the fact that the injected fluid is 
being forced into the formation at a high pressure. The physical property data may comprise an 
actual fracturing pressure, an actual fracturing rate, and an actual fracturing time. 

[0055] The physical property data may be sensed using any suitable technique. For 
example, sensing may occur downhole with real-time data telemetry to the surface, or by delayed 
transfer (e.g., by storage of data downhole, followed by subsequent telemetry to the surface or 
subsequent retrieval of the downhole sensing device, for example). Furthermore, the sensing of 
the physical property data may be performed at any suitable location, including, but not limited to, 
the tubing 1235 or the surface 1224. In general, any sensing technique and equipment suitable for 
detecting the desired physical property data with adequate sensitivity and/or resolution may be 
used. Figure 12 depicts an exemplary embodiment of the present invention wherein the physical 
property data are sensed by a sensing device 1210 resident within well bore 1212. The sensing 
device 1210 may be any sensing device suitable for use in a subterranean well bore. An example 
of a suitable sensing device 1210 is a pressure transducer disclosed in U.S. Patent No. 6,598,481, 
which is hereby incorporated herein for all purposes. In certain exemplary embodiments of the 
present invention, the sensing device 1210 comprises a pressure transducer that is temperature- 
compensated. In one exemplary embodiment of the present invention, the sensing device 1210 is 
lowered into the well bore 1212 and positioned in a downhole environment 1216. In certain 
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exemplary embodiments of the present invention, the sensing device 1210 may be positioned 
below perforations 1230. In certain exemplary embodiments of the present invention, the 
downhole environment 1216 is sealed off by packer 1218, wherein access is controlled with a 
valve 1220. 

[0056] The physical property data is uUimately transmitted to the surface by transmitter 
1205 at a desired time after having been sensed by the sensing device 1210. As noted above, such 
transmission may occur immediately after the physical property data is sensed, or the data may be 
stored and transmitted later. Transmitter 1205 may comprise a wired or wireless connection. In 
one exemplary embodiment of the present invention, the sensing device 1210, in conjunction with 
associated electronics, converts the physical property data to a first electronic signal The first 
electronic signal is transmitted through a wired or wireless connection to signal processor unit 
1222, preferably located above the surface 1224 at wellhead 1226. In certain exemplary 
embodiments of the present invention, the signal processor unit 1222 may be located within a 
surface vehicle (not shown) wherein the fracturing operations are controlled. Signal processor 
unit 1222 may perform mathematical operations on a first electronic signal, fiirther described later 
in this application. In certain exemplary embodiments of the present invention, signal processor 
unit 1222 may be a computer comprising a software program for use in performing mathematical 
operations. An example of a suitable software program is commercially available from The Math 
Works, Inc., of Natick, Massachusetts, under the tradename "MATLAB." In certain exemplary 
embodiments of the present invention, output 1250 from signal processor unit 1222 may be 
plotted on display 1260. 

[0057] , Figure 13 depicts a flow chart of another exemplary embodiment of the methods 
according to the present invention. Li step 1302, the method starts. In step 1304, the method 
creates a fracture layout. Step 1304 may be accomplished using the methods described with 
respect to Figures 3 or 1 1 or any other method for creating a fracture layout. In step 1306, the 
method determines if the fracturing is complete. If the fracturing is complete the method 
proceeds to step 1308, where it ends. If the fracturing is not complete the method proceeds to step 
1310, where a next fracture is induced in the subterranean formation. In step 1312, the method 
receives physical property data from sensing device 1210, the physical property data comprising 
an actual fracturing pressure and may additionally comprise an actual fracturing rate (e.g. a 
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fracturing fluid injection rate), an actual fracturing time, and any surface deformation (eg. a 
fracture-induced surface trough). In step 1314, the method determines if the actual fracturing 
pressure is greater than the predicted fracturing pressure associated with the next actual fracture. 
If the actual fracturing pressure is greater than the predicted fracturing pressure, the method 
proceeds to step 1316, where it modifies the fracture layout and retiuns to step 1306. If the actual 
fracturing pressure is not greater than the predicted fracturing pressure, the method returns to step 
1306. 

[0058] In step 1312, the method according to the present invention receives physical 
property data from sensing device 1210. An example of received physical property data is shown 
in Figure 14, which depicts fracturing pressure (psi) versus time for the creation of six fractures in 
a chert reservoir. The fracturing pressure increases from about 1576 psi for the first fracture to 
about 2600 psi for the sixth fracture. Another example of received physical property data is 
shown in Figure 15, which depicts true vertical depth (TVD) in feet and instantaneous shut-in 
pressure (ISIP) in psi versus measured depth in feet for twenty fractures in a first shale reservoir. 
Yet another example of received physical property data is shown in Figure 16, which depicts true 
vertical depth (TVD) in feet and instantaneous shut-in pressure (ISIP) in psi versus measured 
depth in feet for twelve fractures in a second shale reservoir. 

[0059] Retuming to Figure 13, in step 1316, the method of the present invention modifies 
the fracture layout based on the actual fracturing pressure. In an exemplary embodiment of the 
present invention, the method will reevaluate the fracture layout based on the actual fracturing 
pressure. The method will remodel fractures that have not been induced. The method may use 
the method disclosed in step 304 of Figure 3. The method will substitute the actual fracturing 
pressure for the internal pressure of the next modeled fracture. Based on the reevaluation of the 
fracture layout the method may perform any of the following actions: decrease the number of 
fractures, increase the distance between fractures, or decrease the size of the fractures. For 
example, referring to Figure 2, assume that fracture 206 is the first fracture induced in lateral 104. 
If the actual fracturing pressure associated with fracture 206 is greater than the predicted 
fracturing pressure the method may increase the space between fracture 206 and fracture 204. 
Assuming the actual fracturing pressure is much greater than the predicted fracturing pressure, the 
method may omit fracture 204 entirely, reducing the nimiber of fractures in lateral 104. 
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[0060] The methods disclosed above may be carried out by a computer having a 
processor, a memory, and storage. The methods may be represented as instructions stored in 
software run on the computer. Additionally, the method may be stored in ROM on the computer. 

[0061] Therefore, the present invention is well-adapted to carry out the object and attain 
the ends and advantages mentioned as well as those which are inherent therein. While the 
invention has been depicted, described, and is defined by reference to exemplary embodiments of 
the invention, such a reference does not imply a limitation on the invention, and no such 
limitation is to be inferred. The invention is capable of considerable modification, alternation, 
and equivalents in form and fimction, as will occur to those ordinarily skilled in the pertinent arts 
and having the benefit of this disclosure. The depicted and described embodiments of the 
invention are exemplary only, and are not exhaustive of the scope of the invention. Consequently, 
the invention is intended to be limited only by the spirit and scope of the appended claims, giving 
fiiU cognizance to equivalents in all respects. 
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